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ABSTRACT 

This  report  examines  three  methods  for  maximizing  the  output 
slgnal-to-nolse  ratio  of  a  conventional  passive  detection  system  con¬ 
sisting  of  a  fixed  array,  non-adaptlve  beam  forming,  and  square-law  de¬ 
tection.  Tlie  performance  of  a  non-optlmal  system  is  compared  to  each  of 
the  three  optinum  configurations,  and  the  results  suggest  how  on  unsophis¬ 
ticated,  sub-optimal  system,  which  has  only  marginally  poorer  performance 
than  the  optinum  systems,  can  be  designed. 

The  specification  of  such  a  sub-optimal  system  depends  greatly 
upon  having  intelligent  estimates  of  the  frequency  domain  behavior  of  the 
received  signal;  however,  this  information  appears  only  implicitly  in 
the  solution,  not  explicitly  as  in  the  specification  of  an  Eckart  filter. 
Thus,  one  avoids  the  essentially  insurmountable  practical  problems  of 
atteeqjting  to  build  a  system  matched  to  tlic  signal  characteristics  (tlie 
Eckairt  filter),  while  at  the  same  time  utilizing  the  expected  variation 
in  signal  spectra  to  specify  a  realizable  system  whose  performance  is 
only  slightly  inferior  to  an  Eckart  system. 
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LIST  OF  SYMBOIfi  AM)  ABBREVIATIONS 

It  is  recognized  that  there  will  be  readers  of  this  report 
who  are  not  interested  in  the  mathematical  details,  but  only  in  the 
results  and  their  significance.  To  assist  these  readers,  the  following 
list  of  important  symbols  and  abbreviations  is  included;  references 
are  made  to  pertinent  equations  in  the  text.  This  list,  combined  with 
spot  references  to  the  body  of  the  text,  should  allow  the  interested 
reader  to  obtain  an  understanding  of  the  problem  and  results  by  reading 
sections  1,  2,  and  8. 

w.r.t.  -  "with  resjject  to" 

GNR  -  ubiquitous  throughout  the  text,  denotes  "aignal-to- 
noise  ratio"  and  is  always  written  on  a  power,  not  voltage,  basis  at 
the  point  of  definition  or  measurement;  consequently,  dbs  are  obtained 
by  taking  10  log  of  this  quantity. 

SNRo  -  used  in  two  senses:  (l)  in  a  general  cense,  it  denotec 
the  SNR  out  of  a  square-law  detector  averager  (Eqns.  9  £uid  lO)  and  (2) 
in  numerical  results,  it  denotes  the-  output  SNR  of  a  conventional  bond- 
limited  pre-whitening  syr.tem  (Eqn.  lo). 

SNRo'  (*)  -  system  output  GNR  for  true  Eckart  filter  (Eqns. 

26  and  29):  the  prime  denotec  Eckart,  denotes  tlie  infinite  bandwidth 
associated  with  true  spectral  matching. 

GNRo^  -  cyctem  output  SNR  for  bnndlimiteu  Gckart  filler. 

(Eqn.  3k). 

GNRo'  (§0  )  -  system  output  GNR  for  mismatched  Eckart  filter. 
(Eqn.  U3):  5o  is  the  matching  parameter  (see  helou). 

Si(R,f),  Gi(f)  -  donotes  tlie  signal  power  cpectium  at  the 
input  to  the  shaping  filter.  (Eqn.  1,  Fig.  2)  Because  of  its  depend¬ 
ence  upon  propagation  lose,  the  range  variable,  R,  is  initially  Includ¬ 
ed,  bjt  suppressed  after  section  2  since  the  results  presented  here  do 
not  aepend  upon  it. 
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Nt(f)  -  noise  power  spectrum  at  the  input  to  the  shaping 
filter  (Eqn.  2,  Fig.  2). 

S(R,f),  n(f)  -  signal  power  spectrum  at  the  input  of  the 
square-law  device  (Eqn.  1^,  Fig.  2) 

I  r(f)|*  -  squ!  ire  of  the  trajisfer  function  of  the  shaping 
filter  (Eqn.  13,  Fig.  2). 

A  -  this  symbol  over  a  ixiwer  spectrum,  e.g.  ^(f),  denotes 
single  sided  spectra  (Eqn.  lO).  f.e’ierolly,  double  sided  .spectra  are 
used  in  derivations  while  single  sided  sjiectni  ai-e  used  to  obtain 
numerical  results.  See  Appendix  A. 

Q  -  Power  of  t)ie  frequency  dependence  of  th**  signal  spectra, 

I  .e.,  S'i(f)-Af°'  (Eqn.  12). 

H(f)  -  noise  power  cp'’Ctrun  at  the  input  of  tlie  square-law 
device  (Eqn.  F!g.  .  )• 

3  -  power  of  tlie  frequency  dejiendence  of  tlie  noise  spectra, 
l.e.,  N,(f)«Pf®  (iijn.  11'). 

*  %  S-a  -  referred  to  on  tlie  relative  signol  spectral  slope; 
the  signol  spectra  has  a  slope  of  -10  log  db'oct.  w.r.t.  the  input 
noire  spectra  (Eqn.  I'  ). 

le  -  referr  t  '  tlio  mismatched  Eckart  i’ilter:  5b  is  the 
signal  slope  '*eclgncd  for,  -lilch  is  not  nccecsarlly  the  actual  signal 
rloj'e-  i  (Eqn.  imd  1*2). 

T  -  the  lnt<\Tat  lo:,  time  'iG;u>claVe,l  will.  Ur'  outpu’  am  i>tlilng 
filto;-  (Ecp. .  ‘  and  preceding  discusri  -ri) 

J'l  -  lowf’r  frequency  lltnii  ol'  op<>r’-it  i r.g  l  and:  spectra  are  all 
a;  .u.med  to  he  zero  lelow  thlr  value  (Kqit.  12). 

W  -  the  bar.dwidtl  i  f  Farid  I  Iml led  syrOemr-:  thus,  the  oporalli.g 
land  extends  fi  to  I’j  +  W  (Eqn.  U). 

s  -  rrr.ailzed  lar.dwldti  variable  ( Eqi  .  1") 
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^0  "  the  vuluo  or  bandwidth,  which  maxlmL^eo  tlio  output 
I'.rib  of  n  I’onventionui  pro-whitening;  bandiimitod  oystcm  f  r  each  nignal 
!io|v  1.  (Eqn  iS  and  Fig.  h), 

■Y;(u)  -  i’ui.ctlon  lu’oportionai  to  the  output  Ffili  of  n  conven¬ 
tional  I'andlimited  prc-whitenlng  syctcm  of  bandwidth  u  and  witi  Input 
ifintlve  cignal  rpoctral  rb  jk-  if  t.  "lie  quantity  V^(d)  Y  (ho)  is  the 
ratio  of  tiio  output  oNR  of  a  pir* -whitening  cyatom  of  bandwidth  u  to  the 
outnut  r.'fi'  („f  a  pro-wliltenlng  cyctcn  with  optimum  batidwidth  ho  (Eqns. 
b  and  1  .  Fig.  . 

* 

t  -  multiplicative  factor  dependent  upon  t  al  'i.c,  wliicli  ic 
I’illoaL'.Vf  jf  tiie  pjcrl'omance  of  a  bandlimltcd  pi-e-whi tei.l ng  cyntem 
n  tru<  Fekart  i'llter  (Eqns.  jO  and  ti). 

:'f(u^  -  a  ttinctlon  wiilcii  ic  the  ratio  of  the  output  S^fR  of  a 
band  limited  L’ckart  I'llter  of  boniwldtli  h  to  tlit  output  :’NR  of  Q  true 
Eckart  filter  wl’.h  input  Ignal  c{>oetnii  clop-c  .  ( Bqni  .  3**>  3t>»  'and  <f.). 
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I .  PREUMINAFY  DISCUCGION 


To  oome,  the  optimization  of  a  conventional  cystem  may  seem 
paradoxical  since  tlic  word  conventional  tends  to  imply  non-optlroal  pro¬ 
cessing;  but  the  concept  of  optimi-zatlon  is  multi-faceted,  and  the  term 
remains  ambiguous  until  cystem  configuration  constraints  and  criteria  of 
optimality  are  specified.  In  tills  paper,  an  optimum  system  in  one  which 
maximizes  the  system  output  s Ignal-to-noise  ratio  (f^NR):  from  the  theo¬ 
retical  viewpoint  of  statistical  decision  theory,  this  criterion  is  very 
limited,  and  a  more  satlGfuctoi7  approach  would  be  to  utilize  the  Hayes 
decision  rule  (l,  2,  7,  4,);  however,  most  engineers  arc  more  familiar 
with  the  concept  of  S?fR  than  wltli  Hayes  decision  rules,  and  since,  for 
the  analytical  assumptions  imposed  upon  the  following  work,  the  two  cri¬ 
teria  lead  to  the  same  optimized  configurations,*  the  former  method  is 
chosen . 

Furtlier,  Optimization  is  not  an  end  unto  itself.  Cince  optimum 
cystoms  can  never  be  fully  realized  in  practice,  one  of  the  major  purposes 
of  mat licmat leal  optimization  is  to  determine  the  limit  on  achievable  per¬ 
formance  under  realistic  constraints.  In  view  of  this,  the  primary 
emphasis  of  this  paper  is  placed  upon  the  evaluation  of  a  sub-optimum 
system,  l.e.,  the  relative  cost  (in  dH's  of  output  HfiR)  for  utilizing  a 
non-optinun  cystem  in  place  of  an  optimum  system. 

The  fundamenti.l  constraint  imposed  upon  the  following  work  is 
that  the  cystem  configuration  remains  os  shown  in  Figure  1  for  all  de¬ 
grees  of  optimization.  Tlie  system  consistc  of  an  array  of  lyrdrophonec, 
a  simple  delay  and  cum  beomfonner  (with  constant  multipllvTatlve  spatial 
shading  cocfficlentc  'illowuble),  and  a  signal  processing  unit  consisting 
of  a  square-law  detector,  f.ince  such  a  configuration  forms  the  basic  for 
most  contemporary  systems,  it  will  be  denoted  as  i  conventional  system. 

In  contrast  to  this  conventional  system  are  the  various  "adaptive"  systenr 
wlilch  continuously  modify  hydroplionc  data  on  the  lacic  of  pact  infon-iti- 
tlor.  (5,  o,  ’t  Q);  sucIj  systems  are  still  in  the  conceptual  state,  and 
an.  excluded  from  tliic  analysis. 

•fee  Appendix  D  (l)  for  junpllf ication  of  tills  point. 


1 


NA^T  TR-lUUO 


T))c  opiimlzution  of  the  conventional  cyctem  In  P'lgure  1  io 
further  reotrlctod  to  the  optimal  specification  of  the  shaping  filter, 
lar)!"*  ;  tlie  rirruy  configuration  and  beaffJ*onner  design  are  assumed  fixed, 
and  square-lav  detection  Is  used  throughout.  In  tlie  non-optlmal  cyctem 
which  Is  evaluated,  the  chuplng  filter  Ic  used  only  to  whiten  the  noise 
and  band  limit  the  Inpit  to  the  squnre-lnw  detector;  the  performance  of 
this  system  lo  compared  to  throe  different  optlmur.-  cyntonu  ,  tlie  three 
optlmlzatlonc  rcsnltlnr  ;Yom  different  sets  of  Imposed  contralnts. 

In  the  flrct  optimization.  It  Ic  aceumed  that  tlie  designer  has 
no  fi  priori  knowledge  of  the  slgnol  spectral  chaructcrlctlcs,  except, 
perhape,  a  vigue  fciowledge  of  the  expected  extremo  Unite.  In  this  cace. 
It  Ic  generally  assumed  that  the  best  one  cfm  do  Is  to  wdilten  the  noioe 
and  upproprlatelj  bandllnlt  the  Input;**  thus,  tlilc  optimal  cystem  differs 
from  the  non-optlnul  cyctem  only  to  the  extent  that  tlie  bandwidth  Ic 
(aomewliat  artlflcully)  optimized.  In  fact,  the  procedure  followed  Is  to 
c'/aluote  the  perfona'incc  of  the  non-optlmal  cyctem  as  a  function  of  band¬ 
width,  and  to  denote  that  the  bandwidth  which  maximizes  tlx'  outixit  fHK 
for  1  particular  (but  unknown  to  tlie  designer)  signal  spectrum  as  optimum. 
The  procedure  nay  appear  strained  and  artificial  until  It  ic  reneniiered 
that  the  main  purpose  is  not  to  determine  u  strictly  mathematical  opti- 
niun  system  but  to  evaluate  the  rel/itlvo  perform'inco  of  the  non-optlmal 
cycten.  The  results  will  chow  the  coct.  In  dbc  of  output  r.lfft,  for 
utilizing  an  Incorrect  bandwidth  for  u  p.'irtlcular  signal  spectrum;  per¬ 
haps  more  Important  phyclcally,  the  recultc  also  chow  how  the  performance 
of  u  fixed  bajidwidth  non-optlnval  cyctem  deviates  from  cptlnum  us  the 
InjXit  cpectrum  varies  between  certain  limitn. 

The  cecond  optimization  makes  use  of  tlx-  well  known  llckax^ 
filter  (o,  8,  <)•  In  this  second  case,  It  Ic  uscumed  that  tl  i  designer 
does  I'.ave  prior  knowledge  of  the  signal  spectrum;  this  knowledge  can  be 
effectively  utlllred  to  d«  tennlno  a  form  of  tin*  shaping  filter  which 
maxlni;.es  the  output  .  '.’r.  fuch  a  filter  sp**c!flcatlon  depends  upon  both 
clgnul  and  noise  spectral  cliaricterlct ics,  and  can  be  viewed  uc  the 

••ior  illccucclon,  rt'e  Apjn'ndlx  D  (  ) 
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passive  equivalent  of  the  matched  filter  for  active  cystenc  (l,  10,  11). 

The  non-optisuin  system  performance,  again  as  u  function  of  bandwidth.  Is 
cos^red  to  an  optimal  system  using  such  an  Eckart  filter. 

The  third  optimization  Is  a  modification  of  the  Eckart  result. 

The  assumed  signal  and  noise  spectra  all  approach  zero  asymptotically  as 
the  frequency  approaches  infinity;  thus,  mathematically,  the  Eckart 
filter  has  Infinite  bandwidth.  Physical  considerations,  os  well  as 
curiosity,  leads  one  to  examine  the  effect  of  terminating  the  high 
frequency  response  while  retaining  the  Eckart  filter  characteristic 
over  the  finite  pass-band.  Hie  comparison  made  In  this  third  case  Is 
thus  between  a  bondllmlted  Eckart  filter  and  the  bandllmlted  non-optlmal 
conventional  system,  where  both  systems  are  restricted  to  having  the 
same  bandwidth. 

The  three  cases  illustrate  the  Influence  of  various  filter 
parameters  upon  performance.  The  first  case  demonstrates  the  effect  of 
bandwidth;  the  second,  what  Is  to  be  gained  from  having  and  utilizing 
signal  Information  to  Its  fullest  extent;  the  third  case  coei>lne8  the 
first  two,  and  demonstrates  the  Influence  of  signal  knowledge  In  a  system 
^Ich  similarly  must  be  bandllmlted. 

Thus,  the  results  provide  a  quantitative  estimate  of  the  rela¬ 
tive  value  of  signal  spectral  knowledge;  i.e.,  how  nuch  performance  can 
be  lnq)roved  by  the  Insertion  of  an  Eckoirt  filter  in  place  of  a  simple 
bondllmltlng  pre-whltcner. 

In  order  for  one  to  have  a  keener  appreciation  for  these  results,  the 
next  section  briefly  outlines  the  tactical  passive  sonar  situation  and 
those  factors  which  Influence  the  performance  of  a  conventional  detec¬ 
tion  system. 
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II.  THE  PASSIVE  30NAR  niTUAnON 


For  the  purposes  of  this  trwlysis,  the  passive  sonar  situation 
can  be  depicted  as  shown  in  Kit'.ure  2.  The  fpmntiticc  on  the  left  are 
fenerally  beyond  the  control  of  the  sonar  engineer,  and  nay  be  considered 
as  Inputs.  The  target  which  Is  sought  radiates  noise  of  spectral  con¬ 
tent  T(f).  The  signal  spectrum  received  at  the  detecting  platform  is 
T(f)  as  modified  hy  Lh'*  propagation  loss  H(F,f).  Masking  the  target 
radiated  noise  are  both  ••nvl ronmental  and  platform  generated  (self)  noise; 
it  .'.s  assumed  that  tnes**  two  noise  processes  ore  independent  so  that  their 
p)owsr  spectra  combin»-  art  !  i  t  ivt- ly.  This  total  noise  is  reduced  by  the 
array  gain,  so  that  the  resultant  signal  and  noise  spectra  into  the  shap¬ 
ing  filter  are* 


S.  Cf  f)  -- 
N,(f)  = 


TCf)  H(R,4) 

— Af^i 


tl) 

(2) 


The  spectral  characteristics  of  the  quantities  on  the  right  of 
equations  (1)  and  (2)  can,  of  course,  be  rather  complex;  however,  in 
regiona  of  interest,  the  following  trends  are  generally  exhibited:** 


T(f):  spectral  slope  of  approx iauitely  -5db/oct. 

H(R,f):  highly  variable  spectral  slope,  depending  upon  mode  of 
propagation,  range,  and  e.nvironmental  conditions;  typical  slopes  vary 
between  Odb/oct  and  -20dl'/oct. 

Ns(f),  spectral  slope  of  approximately  -Jdb/oct. 

A(f);  spectral  slope  of  approximately  +bdb/oct  for  planar  and 
volumetric  arrays,  +3db/oct  for  llae  arrays*** 

*This  formulation  as^tiines  the  target  is  on  the  ma.lor  response  axis  of  the 
baam;  an  off  axis  target  would  naturally  have  its  spectrum  modified  by 
the  frequency  dependouoo  r,r  tiic  beam  response. 


••These  values  are  presented  only  to  give  a  physical  feeling  for  the 
problem  being  treated.  They  can  be  ignored  from  a  mathematical  standpoint, 
since  the  problem  can  be  viewed  simply  as  an  optimization  exercise.  That 
the  stated  trends  are  representative  can  be  deterraiued  ly  examining  refer¬ 
ences  12  and  13 . 


••♦Arrays  assumed  to  be  in  their  aperture  response  region,  i.e.,  less 
than  half  wave  length  spacing  of  ’he  elenents. 


f 

) 

1 
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Figlire  2  Passive  Sonar  Situation 
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Ao  Implied  from  the  preceding,  the  spectral  characteristics 
of  :'j(R,f)  exhibit  n  lilglier  degree  of  variability  than  those  of  Nj(f) 
by  virtue  of  the  proijugatlon  loss,  H(R,f).  Consequently,  it  is  only 
natural  to  examine  what  effect  this  variation  has  upon  the  performance 
of  a  passive  detection  system. 
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III.  DfcyiVATION  OF  CONVEJmONAL  JYSTEM  K^UATION 

In  order  to  proceed  with  the  performance  evaluation,  it  is  nec¬ 
essary  to  determine  the  relationship  between  the  system  output  SNR  and 
the  signal  and  noise  spectra  into  the  square-law  device.  Compared  to  the 
remainder  of  this  paper,  this  section  is  rather  mathematical.  The  reader 
stiould  have  a  basic  understand iriK  of  Fourier  transfoms  and  their  manipu¬ 
lation  ns  well  as  a  familiarity  with  the  fundamental  concepts  of  si(^nal 
processinr;  the  second  cliaptcr  of  Wi.K>dwnrd 's  monograph  provides  an  adequate 
coverage  of  this  ma’erial  (lO).  Tliose  not  so  endowed  can  skip  to  the  next 
section  with  11* 'le  loss  iti  continuity. 

This  section  is  concerned  with  only  two  I  locks  ol'  the  processinr 
functions  shown  In  Flt'ur*-  2;  these  are  the  square-law  device  and  the  out¬ 
put  smoothing.  filt'?r,  shown  in  Firure  3-  The  input,  consists  of 

either  sirnal  plus  noise  or  noise  alone: 

Z  \ 

I  ^(t)  [it) 

Assume  ttmt  s(t)  and  n(t)  ure  members  of  a  zero  mean  stationary  Gaussian 
random  process.  Then  the  correlation  function  at  the  output  of  the  square- 
law  device  is  (lU,  p.  255). 


Further  assume  th.ut  s(t)  is  statistically  independeri*  of  n(t),  and  that 
they  have  variances  Sj*'  and  then 

\  O'- 

L  «•«  *  .•>  if.’-C*') 

Realizinir  that  the  power  spec’ral  density  of  a  random  process  is  the 
Fourier  transform  of  its  correlation  function,  and  that  multiplication 
in  the  t  domain  corresponds  to  convolution  in  ’he  frequency  domain,  the 
fxjwer  spectral  density  of  y(t)  is 

Vs'  +  /vcvj 

<  or 


(Hi) 

(Mfc) 


(5a; 


£(f )  +7  N(i)  i)  no 


(£{>) 


I 

J  I 

N 

J 
I 
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where  ^(f)  Is  the  Dirac  delta  functlODf  S(f)  ajid  N(f )  ore  the  signal 
signal  and  noise  power  spectra  into  the  square-law  device,*  and  the  syidbol 
■  denotes  convolution. 

The  wave  form  y(t)  is  averaged  by  a  filter  with  transfer  function 
H(f );  thus  the  output  power  spectrum  is 

Sz(f)  -  Sy(f)  |H(f)|*  U) 

The  results  of  equations  (3)  and  (6)  contain  the  information 
needed  to  determine  the  system  output  SIfR,  which  must  now  be  defined.  The 
output  signal  of  interest  is  the  shift  in  the  average  (d-c)  voltage 
level  of  the  output  due  to  the  presence  of  signal.  Assuming  that  the 
Input  signal  and  noise  spectra  do  not  contain  In;>ul8e8  (l.e.,  ignoring 
the  effect  of  lines),  and  assuming  |h(0)|*  1,  the  output  d-c  voltage  %d.th 
signal  present  is  Oj  +  ON*(the  6(f)term  from  equation  5a).  Since  the  d-c 
level  with  noise  alone  is  (equation  5b),  the  output  signal  voltage  of 
interest  is  Oj  .  Returning  to  a  power  basis  in  order  to  write  SNR's, 
the  output  signal  power  lo  a j  ^ . 

The  continuous  portion  of  SzCf)  is  the  "ripple",  or  noise, 
present  at  the  smoothing  filter  output  which  tends  to  obscure  the  shift  in 
d-c  level  due  to  the  pwsence  of  signal.  There  are  two  possible  definitions 
of  output  noise:  one  would  Involve  using  the  continuum  portion  of  the 
spectrum  in  equation  5<^>  l.e.,  the  ripple  due  to  the  presence  of  both 
signal  and  noise  at  the  input;  the  other  Involves  the  continuum  portion  of 
the  spectrum  in  equation  5b,  l.e.,  the  ripple  due  to  the  presence  of  noise 
alone  at  the  input.  For  a  complete  description  of  processor  performance, 
both  quantities  ore  needed;  however,  in  most  passive  sonar  situations  of 
interest,  the  input  SNR«1,  and,  in  this  case,  the  presence  of  a  signal 
has  a  negligible  effect  upon  the  output  variation.  Because  of  this  effect, 
as  well  as  the  resulting  mathematical  simplicity  and  the  degree  of  arbi¬ 
trariness  associated  with  any  definition  of  SNR,  the  output  noise  is 

*The  range  variable,  R,  has  been  suppressed  since  it  exerts  no  Influence 
upon  the  following  results. 


) 


J 
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defined  to  be  the  output  verletlon  due  to  nolee  alone  at  the  Input.* 

Thus,  the  output  SNR  is,  from  equations  5  and  6, 

s s  »»**  {  2  f  ^(0  a I  j in) 

This  expression  can  be  siispllfied  a  great  deal  by  one  further 
aasuarption  often  satisfied  in  practice.  Assume  that  the  output  filter 
bandwidth  is  very  narrow  coBq;>ared  to  the  input  bandwidth  of  the  system; 
under  these  conditions,  the  auto-convolution  of  N(f)  will  have  a  very  broad 
peak  in  the  region  of  f-0  %rhen  compared  to  ]H(f)p.  The  value  of  this  peak 
will  be  the  auto-convolution  evaluated  at  f>0  i.e., 

m 


since  N(-S)  «  N($),  i.e.,  the  power  spectra  of  real  functions  are  real  and 
even.  Utilizing  this  asstuption  that  the  output  filter  response  is  suffi¬ 
ciently  narrow  so  that  |H(f)|^  0  before  the  value  of  R(f)  sdf(f)  has  varied 

appreciably  from  its  value  at  f>0,  and  further  recognizing  as  the  input 
slgial  power,  equation  7  can  be  tiffproxiMMtti.  by 


SVA.il  _ _ 


Equation  8  is  a  rather  general  expression  for  the  output  SNR  of  a  square-law 
detector- averager.  Die  only  restrictive  assuaptlons  used  to  derive  this 
equation  are 

<^is  definition  of  output  SNR  corresponds  to  the  following  statistical 
definition: 

[r7s+)!^)  ■ 


where  Z  is  the  output  mesui,  is  the  output  varleuice,  and  N  and  S  +  N 
denote  the  presence  of  noise  alone  or  signal  plus  noise,  respectively. 
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(1)  tb«  iapat  ii^aal  and  nolM  ar*  Mabara  of  atatlooanr  OauaaUn 
raadoa  proeetaaa,  indtpendant  of  ona  another,  and 

(2)  tha  Ixq^t  nolaa  la  *Snda  band"  eoqparad  to  tha  output  flltar. 

For  the  remainder  of  thia  paper,  the  output  arooothlng  filter 
la  ehoaen  to  be  a  perfect,  finite  time  averager  with  Inpulae  reaponae 


0<t<T 

( \stMjka^c 


Than  (10,  13) 


Thua,  for  thia  output  filter,  the  expreaaion  for  the  output  SNR  may  be 
written  aa 


IZ 


If  aingle  aided  power  apectra,  S(f)  and  N(f),  are  used  in  place 
of  tha  double  aided  S(f)  and  lf(f),  the  reatHting  expreaaion  ia 


S/VR.  •=  T 


aa  ahovn  in  Appendix  A. 

If  the  nolae  la  white  over  a  band  W,  it  la  a  straightforward 
proceaa  to  ataofw  that  equation  10  (or  9)  reduces  to 

SNR,=  TW(A;^  (III 

where  S  and  N  are  the  signal  and  noise  power  (total  power,  not  spectral  levels) 
Into  the  square-law  device. 

*8ee  Appendix  B  ^or  a  generalization  of  equation  (U). 
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IV,  CASE  1;  BAWEWIETH  OPTIMIZATION  IN  WHITE  NOISE 


In  order  to  obtain  numerical  results,  some  functional  form 
must  be  assumed  for  the  Input  signal  and  noise  spectra,  S|(f)  and  Ni(f) 

In  Figure  2.  Since  many  of  the  quantities  in  the  sonar  equation  can  be 
conveniently  represented  by  a  db/octave  slope  over  the  region  of  interest, 
the  following  dependence  is  assumed: 


f,  f  VV 

t,  f  “ Ir*/ 

ujk9 


(J2‘) 


The  spectra  are  thus  assumed  to  be  shai*ply  bandlimlted,  and  to  be  charac¬ 
terised  by  constant  db/octave  slopes  w'thin  the  band  (f:,  fi  +  W).  This 
Is  an  oversimplified  representation,  particularly  for  Si(f)  because  of 
the  propagation  effects  mentioned  earlier;  however,  when  tempered  by 
reason,  the  results  from  such  spectra  shed  light  upon  passive  system 
performance . 


The  spectra  Into  the  square-law  device  are 

To  whiten. the  noise,  the  shaping  filter  must  be  of  the  form* 


so  that 


A-f 

P(f)  =  & 


n 


OJo) 

(JH) 

(/S’i) 


where  5=0-0  is  the  input  signal  spectral  slope  with  respect  to  the  in¬ 
put  noise  spectral  slope,  l.e.,  the  slope  of  the  signal  spectrum  Into 
the  square-law  device. 

The  direct  application  of  equation  (lO)  to  the  spectra  in 
equation  (15)  results  in  , .  .,.Z  ^  ) 


rfir 


,  S  =  1 


(lio) 

t'4*) 


*The  arbitrary  multiplicative  constant  is  chosen  to  be  unity  for 
simplicity. 
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where 


(/•)  =  ;!:  If' .  !♦  I 
(a)  -  .?*' 


(•7.j 

07i) 


1 

1 


and  (Jk  w/fi  is  the  normalized  bandwidth  variable,  which  gives  the  sys¬ 
tem  bandwidth  as  a  multiple  of  fi.  Thus,  u  =  1  corresponds  to  an  octave 
jand. 

In  any  given  sonar  situation,  wltii  an  assumed  constant  rela¬ 
tive  signal  spectral  slope  T,  A,  B,  fi ,  and  ^  are  all  constant.  Thus, 
the  describe  the  relative  SNR  performance  as  a  function  of  band¬ 

width.  From  the  expressions  In  equation  17,  It  can  be  seen  that  the 
possibility  exists  for  V^(k*)  to  have  a  relative  maximum  at  some  value 
u  =  depending  upon  the  value  of  5.  Considering  the  5  /  1  case,  set¬ 
ting  0  results  in  the  following  equation  for  determinant  mq: 

n-/.(2«-i) .  (/y-A)'  ,  5#i  00 

Thus  if§  =  2,  M  =  1;  consequently,  for  a  relative  signal  spectral 
slope  of  -6db/oct,  the  bandwidth  which  msocimlzes  the  output  SNR  is  an 
octave  wide. 

there  does  not  always  exist  an  optimum  finite  bandwidth,  if 
&  »  0,  i.e.,  the  signal  spectrum  into  the  square-law  device  is  also 
white,  then  equation  l8  results  in  Mo  »  0,  which  is  obviously  not  the 
relative  maximum  sou^t.  Referring  to  equation  (11),  if  both  signal  and 
noise  are  white,  then  the  SNR  into  the  square-law  device  Bqn.  ll) 

is  a  constant  as  W  is  varied,  and  equation  (ll)  Indicates  that  the  out¬ 
put  SNR  will  Increase  lineeurly  with  W.  Thus,  the  bandwidth  which  would 
maximize  the  output  SNR  in  this  case,  ^  =  0,  would  be  infinitely  wide. 
Naturally,  such  a  spectrum  will  not  occur  physically  as  various  limiting 
mechanisms  will  eventually  cause  the  high  frequency  response  to  fall 
off;  however,  mathematically,  no  finite  uo  exists. 

Looking  at  equations  (11)  and  (l8),  it  can  be  seen  that  an 
infinite  bandwidth  solution  will  result  until  some  critical  value  of 
is  used  for  the  relative  signal  spectral  sloiw.  Thus,  in 
range  0  ^  ^  ^  io ,  even  though  the  signal  spectrum  is  falling  off 
compcu^d  to  the  noise  spec  trim,  and,  as  a  consequence,  the  SNF  into  the 
square-law  device  is  decreasing  as  W  is  increased,  signal  processing 


Ih 
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(the  TW  portlu.-.  of  Eqn.  11)  nore  than  cornpensatea  for  the  decreaelng 
Input  SIJR.  Intuitively,  one  would  suspect  the  value  of  ?o  to  be  ^ 
(corresponilng  to  a  relative  signal  spectral  slope  of  -l,5db/oct);  for 
once,  intuition  proves  correct.  For  ?  =  equation  iB  shows  that  no 
relative  maxitiiur  exists,  an'  *te  infinite  bandwidth  solution  again  re¬ 
sults.  Appendix  C  demons trates  that  5  j>  Is  the  critical  value,  and 
that  for  values  of  ?  Incremental!;,  larger  than  finlt«!  bandwidth 
solutions  result. 

One  -ouli  use  equation  iB  to  letermlne  optimum  bandwidths  for 
various  valuer  ol  ?(!  /  1,  i  >  ‘);  however,  it  is  not  only  the  optimum 
value  of  bandwidth  i,hat  is  of  concern,  tut  also  the  penalty  that  is  paid 
for  not  ojh-rating  at  the  optimum.  This  is  simply  the  evaluation  of  the 
non-opt ImuTT.  system  discussed  in  Section  1.  Such  an  evaluation  can  be 
accomplished  by  examining  a  normalized  as  a  function  of  ^  for 

various  values  of  ^  In  the  rci^ion  where  optimum  Vandwidths  exist.  This 
normalized  Vr(-),  given  h;,  Y.( >  Is  plotted  in  Figure  U  for  5 

••  s 

equal  to  1  througli  o  In  integer  values. 

Referring  to  equation  10,  it  can  be  seen  that  Yr(^^)/Y^(^*o) 
is  simply  the  ratio  of  the  output  SNR  at  a  normalized  bandwidth  to 
the  output  SNR  for  the  optimum  bandwidth,  >4o,  corresponding  to  the 
particular  value  of  1.  The  quantity  Y^(^)/Y|{i*o )  can  thus  be  conven¬ 
iently  referred  to  as  the  relative  output  SNR  since  it  measures  the  SNR 
out  of  a  non-optimum  oystem  relative  to  the  output  SNR  that  would  be  ob¬ 
tained  had  an  optimum  bandwidth  been  employed. 

The  graphs  in  Figure  b  are  intuitively  satisfying  since  the 
peaks  are  rather  broad,  indicating  little  loss  associated  with  picking 
a  non-optimum  band.  Often,  an  octave  band  l)  is  used  lii  acoustic 

analysis  systems.  Figure  ^  shows  that  for  sit-n^al  spectra  into  the 
souare-law  device  with  slopes  between  -3db/oct  and  -iBdb/oct,  the  loss 
in  output  SNR  associate!  wi‘h  an  octave  land  is  less  than  3db.  Thus, 
Figure  b  would  indicate  that  an  octave  's  a  good  compromise  band.  It 
must  be  stressed,  however,  that  this  section  has  dealt  with  extremely 
simple,  and  consequent!;',  somewhat  artificial,  spectral  characteristics. 
The  results  are  enlightening  and  meaningful,  but  simply  In  pointing  the 
way  toward  a  solution;  they  do  not  provide  a  comprehensive  answer  to  the 
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Figure  i*-  Relative  Output  SNR  for  Bandllmlted  Prevhltening  System 
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passive  sonar  detection  problem.  Besides  the  simplified  spectral 
characteristics,  no  mention  has  been  made  about  determining  f i ,  basic¬ 
ally,  specifying  the  location  of  the  operating  band.  Such  a  specifi¬ 
cation  is  strongly  dependent  upon  the  basic  sonar  inputs  of  Fit;  ire  2  -- 
the  target,  noise,  and  propagation  characteristics  --  as  well  as  physical 
limitations  upon  the  array  size  and  available  equipment  space  for 
beamforming  and  processing.  The  consideration  of  all  these  factors 
would  lead  to  a  total  system  design,  which  is  well  beyond  the  scope  and 
intent  of  t’nis  article;  however,  some  further  general  comments  about 
the  influence  of  these  quantities  upon  system  bandwidth  specification 
will  be  made  in  the  concluding  section. 
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V.  CASE  2:  THE  ECKART  FILTER 


The  previous  section  examined  what  might  be  termed  a  theoreti¬ 
cally  crude  optimization  scheme:  the  shaping  filter  was  constrained 
to  be  a  noise  pre-whitening  filter  only,  and  the  optimization  process 
was  performed  only  on  bandwidth.  In  this  section,  the  optimization  is 
more  sophisticated,  making  use  of  signal  spectral  knowledge  to  deter¬ 
mine  the  best  form  of  shaping  filter.  The  result  is  the  well  Eckart 
filter,  given  by  ^ \ 


|r(f)l 


s.q) 


To  show  that  this  is  the  optimum  filter,  equation  9  will  be 
maximized  by  a  common  procedure  (8):  assume  a  general  shaping  filter 
response  given  by 

r(f)  .=  ro(f)  +  ev(f)  (7.0) 

where  ro(f)  is  the  optimum  value  of  r(f)  sought,  Y(f)  is  an  arbitrary 
frequency  function,  and  e  is  any  scaler.  Then  the  system  output  SNR 
is,  from  equations  9f  13)  and  20, 

■°*  dr 

If  lo  is  truly  the  optimum  filter  response,  then^  should  be  zero  when 

e  is  set  equal  to  zero,  i.e., 

^  {  D,  C22) 

'  K  C ^  r/'m] 

where  Do  and  No  represent  the  denominator  and  numerator  '•f  equation  21 
evaluated  at  e  =  0.  Further  manipulation  results  in 


[  (vr/+  x*n][  s-{f)  -KA//0)lrJ^J  df  =0 


(2  3) 
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where  K  =  .  The  only  way  that  equation  23  can  be  valid  for  all  values 

of  y{f),  which  is  an  arbitrary  function,  is  if 


r  ff)  I  ^ 

K  ^;(f) 


Since  K  is  a  scalar  constant,  it  acts  merely  as  a  scale  factor  which  can 
have  no  effect  upon  the  optimum  specification  of  r(f );  thus,  it  can  be 
arbitrarily  chosen  as  unity,  resulting  in  the  specification  of  equation  I9. 

Using  the  Eckart  filter  for  Jr(f)j*,  the  spectra  into  the  square- 
law  device  become 

[TvZfjJ 


N(^)  - 


V  (f-) 
NiCO 


(2  s-l>) 


By  substitution  into  equation  10  (or  9),  the  output  SNR  takes  the  simple 
form  ^  </.  i  1  ? 

S/v^r-)-  t/  1%^)] 

n 

Utilizing  the  same  spectral  forms  assumed  in  Section  4,  but  ex¬ 
tending  the  high  frequency  response  to  infinity,  the  input  signal  and 


(24) 


noise  spectra  are 


we  have 


>f’fl  ] 

^  o<  <0 

3 

.f>f,  l 

,f^f,  J 

p<o 

(z  nt) 

t 

^f>f, 

(i?) 

mW  ^ 

where  §  =  3-0;,  as  in  Section  4. 

Substituting  equation  28  into  26  and  integrating,  the  output  SNR 
for  an  Infinite  bandwidth  Eckart  filter  is 


s<v^Y‘»)=  (if  (7^;  ^ 
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Referring  to  equation  l6,  the  performance  of  the  non-optiimun 
conventional  system  with  normalized  bandwidth  m>  can  be  compared  to  the 
Eckart  filter  system.  The  relative  loss  in  output  SNR  of  the  non-optimum 
system  is  conveniently  expressed  as 


where  is  a  multiplicative  factor  depending  on  §  alone,  amd  given  by 


Thus,  graphs  of  SNRo /SNRq ^  (“)  exhibit  the  same  bandwidth  dependence  as 
Figure  Uj  the  curve  corresponding  to  each  value  of  |  is  simply  moved  down  by 
a  factor  of  which  is  a  constant  for  each  value  of  The  results 

for  equation  30  are  shown  in  Figure  5*  It  can  be  seen  that  the  additional 
"signal  shaping"  loss,  i.e.,  the  loss  in  addition  to  the  bandwidth  loss 
shown  in  Figure  U,  amounts  to  some  1  to  2  db,  with  the  higher  loss  being 
associated  with  the  shallower  relative  signal  spectral  slope,  §  =  1. 

Figure  5  shows  that  an  octave  wide  pre -whitening  system  loses  between  1 
and  4  db  in  output  SNR  compared  to  an  Eckart  "matched"  system  for  relative 
signal  slopes  between  -3db/oct  and  -l8db/oct. 
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Figure  5  Con^arison  of  Bandlimlted  Prewhltening  System  (SMRj)  to  Eckart  Filter  System  (SNRo(°°)) 
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VI.  CASE  3:  THE  BANDLIMITED  ECKART  FILTER 


The  last  optimization  to  he  examined  arises  basically  from 
curiosity:  one  wonders  just  how  sigiiificant  the  higher  frequencies  are  in 
the  specification  of  the  Eckart  filter.  In  this  section,  the  performance 
of  both  the  infinite  bandwidth  Eckart  filter  and  the  pre -whitening,  band- 
limited  conventional  system  will  be  compared  to  a  bandlimited  Eckart  filter, 
defined  as 

c 


|r(f)|^  - 1 


C'p) 


t  (it 


The  system  output  SNR  for  such  a  filter  follows  from  equation  26: 


S<uR, 


•f 


Jf 


Utilizing  the  spectral  characteristics  from  equation  27, 


where 


Vlj  (/a)  ~  /-  </'*/')  ^ 


(33) 


(3M) 


and  the  logarithmic  solution  for  f  }  is  ignored. 

By  comparisoii  with  equation  29,  it  can  be  seen  that  the  ratio  of 
the  bandlimited  Eckart  filter  to  infinite  bandwidth  Eckart  filter  output 
SNR's  is  simply  to  factor  in  equation  3?^  i.e., 

The  comparison  to  the  pre-wliitening  bandlimited  conventional  system  can 
be  rapidly  obtained  from  the  following  irenipulation: 


(3(^) 


][  C»-)7 


(37> 
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Both  factors  on  the  right  of  equation  37  have  already  been  determined, 
viz.,  from  equatione  30  and  36,  so  that 


S/v«/ 


Since  we  are  only  concerned  with  those  cases  for  which  optimum  bandwidths 
exist  for  ■Y|.((-),  the  value  of  §  can  be  restricted  to  being  greater  than  i, 
and  no  perplexing  mathematical  difficulties  are  encountered  in  the  above 
manipulations . 

Equation  36  ic  plotted  in  Figure  6;  naturally,  T1^(M')>  the  loss 
associated  with  bandlimiting  the  Eckart  filter,  approaches  Odb  as  |i  gets 
large,  the  approach  being  fastei'  as  §  increases.  The  use  of  a  bandlimited 
Eckart  filter  with  p.  =  uq,  the  "optimum"  bandwidth  for  each  §  from  Figure 
4,  results  in  I-  to  1  db  of  loss  compared  to  the  infinite  bandwidth  Eckart; 
recall  that  for  the  conventional  pre-whitening  system  with  n  =  no ,  the  loss 
was  from  1  to  2  db  (Figure  5)* 

Equation  38  is  plotted  in  Figure  7-  Again  the  results  agree  with 
intuition.  Where  both  the  Eckart  and  pre-whitening  filters  have  narrow 
bandwidths,  they  result  in  approximately  the  same  performance.  For  large 
bandwidths,  the  factor  ’’-(n)  in  equation  38  approaches  unity,  and  the 
bandlimited  Eckart  cut- performs  the  pre-whitening  system  by  the  same  amount 
as  the  infinite  bandwidth  Eckart  shown  in  Figure  5  • 
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TABLE  I  Comparison  of  Mismatching  Losses 


(a)  Systems  designed  to  optimize  performance  at  §  =  1 


Actual  value 
of  § 

Loss  of  pre -whitening 
system  w.r.t  optimum 
pre-whitening  system 
(db) 

Loss  of  pre-whitening 
system  w.r.t.  matched 
Eckart  system 
(db) 

Loss  of  mismatched 
Eckart  system  w.r.t. 
matched  Eckart  system 
(db) 

1 

0.0 

-1.9 

0.0 

? 

-2.0 

-3.4 

-1.3 

3 

-t.L 

-5.5 

-2.6 

h 

-6.4 

-7.5 

-3.6 

5 

-7.5 

-8.7 

-4.4 

r 

U 

-8.8 

-9.7 

-5.2 

(b'l  System 

s  designed  to  optimize 

performance  at  §  =  3 

Loss  of  pre-whitening 

Loss  of  pre -whitening 

Loss  of  mismatched 

nf  c 

system  w.r.t.  optimum 

system  w.r.t.  matched 

Eckart  system  w.r.t. 

Ui  ‘3 

pre-whitening  system 

Eckart  system 

matched  Eckart  system 

(db) 

(db) 

(db) 

1 

-2.5 

-4.4 

-2.5 

? 

-0.4 

-1.6 

-0.4 

3 

0.0 

-1.2 

0.0 

U 

-0.2 

-1.3 

-0.2 

5 

-0.6 

VO 

• 

1 

-0.5 

u 

-1.1 

-2.1 

-0.7 

IF 
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V'JII.  nONCLUDINC  DIGCUGGION 

Tlie  specific  I'esuLf.'  of  this  report  are  contained  in  Figures  4 
tlirougli  8.  Althougli  ;',omc  explanation  and  interpretation  of  these  graphs 
.-.'as  included  in  cacli  of  the  indi\/idual  sections  where  they  initially 
appearoo,  this  section  v,'ill  expand  and  clarify  the  results  further.  For 
convenience  In  reference,  the  five  figures  are  reproduced  at  the  end  of  this 
section . 

Figure  4  sl.ows  that  the  f)enaity  for  using  ;ni  incorrect  bandwidth 
is  n  't  severe,  tl'.ough  an  intelligent  choice  of  bandwidth  can  reduce  tliis 
loss  by  a  significant  amount.  As  an  example,  assume  that  the  signal  spec¬ 
tra  are  expiected  to  Ivvo  relative  slopes  varying  between  -3  and  -l8db/oct. 

If  tlie  system  used  a  bandvidtli  corresponding  to  a  -  8,  the  loss  (w.r.t. 
optimum  u's)  in  outijut  f-NR  would  vary  between  0  and  lOdb ,  If  the  system 
were  redesigned  so  tliat  o.  =  0.78  used,  the  loss  would  be  betweeri  0  and 

db - an  8  db  improvement  in  worst  case  system  performance  achieved  very 

simply  by  choosing  the  "correct"  bandwidth.  If  the  occurrence  of  the  various 
i'ignal  slopes  wore  not  equi -probable,  tlie  bandwidth  choice  might  be  altered. 
For  example,  if  slopes  of  - and  -odb/oct  occurred  90^  of  the  time,  then 
one  might  choose  u  ''w  2  since  the  "majority"  .signals  v/oalri  suffer  less  than 
^  db  of  loss.  With  a  u  of  2,  the  v/orst  case  loss  is  about  5.^'  db  more 
tlian  tlie  worst  case  loss  associated  with  u  =  0.75'  hut  heie,  the  worst 
cast?  loss  is  sufi’ered  by  the  least  probable  signals,  while  with  M'  = 
the  worst  case  loss  of  2  n.b  would  be  .uffered  by  one  of  tlie  more  probable 
signals.  All  in  all,  of  coux'ce,  the  final  choice  is  a  iiiattex’  of  the 
designer's  .ludgement — no  "rules"  can  be  given;  hov/over.  Figure  4  does 
contain  information  which  should  bo  useful  in  making  such  a  Judgement. 

figure  5  sliows  tie  ndditi  nial  losses  entailed  by  not  fully 
ut ilising  spectral  information  ab  ut  the  signal.  Notice  that  if  the 
region  of  concerii  were  again  1  ^  r  n  pjje  niinlmum  worst  case  loss  i'or  a 
conventional  bandlimitod  pro-whitening  system  no  longer  occurs  at  m-  =  0.78> 
but  at  K  --  t'.9;  tliG  change  in  the  criterion  of  comparison  has  affected 
the  choice  of  "correct"  bandwidth,  although  not  by  a  significant  araout.t. 


30 


NAFT  TR -141+0 


Naturally,  as  tlie  regiJii  of  signal  interest  changes,  so  does  the  choice  of 
bandwidth;  if  the  possible  signal  slopes  lie  between  j  ^  §^5  (between 
and  -I5  db/oct),  the  bandwidth  resulting  in  the  smallest  worst  case  loss 
would  be  0.4,-;.  For  the  range  of  signal  slopes  pi-  'ted,  the  loss  over 
and  beyond  the  bandwidth  loss  of  Figure  4  amounts  to  some  1  to  2  db . 

I'igures  0  and  '(  sliow  how  the  swo  preceding  systems  compare  to 
a  bandlimited  Eckart  fiitrT.  l''igure  6  indicates  that  there  is  little  use 
to  extend  the  bandwiJtl:  ol‘  an  Eckart  filter  beyond  about  2  ,  where  the 

Po  are  opt.imiim  bandwidths  determined  from  Figure  4.  Figure  7 
i.ndicntes  tlia.t  a  nari'ow  bandwidth  Eckart  filter,  witli  a  po/.',  is  of  little 
value  since  its  per''ox’'.iance  is  nearly  iJeiatical  to  the  simple  bandlimited 
pre-wliitening  system. 

Figure  3  sho-s  the  losses  Involved  in  using  .a  mismatclied  Eckart 
filter,  i.e.,  using  a  filler  designed  for  to  when  tlie  signal  spectrum 
actually  has  a  slope  t-  I'y  comparison  wltli  Figures  4  and  '■j,  it  can  be  seen 
that  the  Eckart  i'ilter  gTonerully  exliibits  less  mismatching  loss  than  the 
pre-whitening  bandlimited  system.  If  the  designer  were  constrained  to  a 
single  shaping  i'ilter  specification,  i.e.,  he  could  not  implement  a 
matching  (Eckart)  filter  for  each  and  every  signal  spectra,  tlien  a  mis¬ 
matched  Eckart  filter  with  to  *'«'  2.2  would  minimize  tlje  worst  case  loss  for 
1  -  i  ^  o.  The  performance  of  such  a  mismatched  Eckart  system  is  compared 
to  the  minimum  worst-case,  bandlimited,  pre-wiiitening  system  in  Table  II. 
Wlien  the  performance  of  both  systems  is  compared  to  a  m.atched  Eckart  filter, 
the  mismatched  Eckart  is  superior  by  some  1  to  2  db . 

Tlius,  the  following  conclucion.s  can  be  drawn  from  the  numerical 
results  contained  in  this  paper.  For  signal  spectra  which  have  a  constant 
slope  w.r.t.  noise  "'f  -lU  log  2"  db/oct  above  some  lower  limit  of  fi  cps: 

(1)  Any  Eckart  filter  implementation  need  not  have  a  bandwidth 
greater  than  2  so  .  where  Hie  Po  are  def.ermined  from  Fi  aire  4. 

(2)  An  Eckart  filter,  for  any  particular  value  of  i,  improves 
performance  by  1  to  I  db  over  an  optimized  bandlimited  pre-whiloning  sys- 
t.em;  aith  ugh  noP  totally  insignificant,  tlie  improvement  Is  marginal. 

i  ■])  h  mismatclied  Eckart  filter  exhibits  less  loss  than  a  mis- 
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TABLE  II  Comparison  of  Minimum  Worst  Case  Design  to  Optimum  Systems 


1 

fiOor;  of  mitilnaii!;  v/orol 
case  pre-wh  !  ten!  iifj 
r.ystem  v.r.t.  opti'mn:, 
pre-wh i  ten  1  nr 
(el) 

IjOse  of  minimum  wor.oL 
case  pre-wh iteninr 
system  v.r.t,  tnatcheh 
Eckart  syslett; 

(ai) 

r,oss  of  minimum  worst 
case  mismatched 

Eckert  system  v.r.t. 
matclied  Eckart  r;vstem 

(lb) 

1 

-1.88 

-3.38 

-1  .If' 

n 

< 

o 

o 

1 

-]  /'('■> 

-0.03 

i 

-0 . 10 

-1.37 

-0.18 

-0 . 6l 

1 

b 

O 

-O.oi 

3 

-1.2u 

-r.oB 

-l.l7 

'J 

-1.88 

-3.38 

-1.1+9 
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maicVied  pre-whiteninfi;  system.  If  a  minimum  worst  case  loss  criterion  is 
used  to  specify  p,,  the  bandwidth  of  the  pre-whiteninr;  system,  and  5o  j  the 
"matchinf^:  parameter"  of  the  Eckart  filter,  the  Eckart  desicin  is  superior 
by  1  to  2  db.  Au:ain,  the  improvement  appears  marfunal. 

The  preceding'  results  and  conclusions  are  strictly  valid  only 
for  the  assumed  simple  spectral  characteristics;  iiov/ever,  the  overall  trends 
exhibited  lead  one  t,o  formulate  hypothesized  c;eneralizations .  The  impli¬ 
cations  for  conventional  system  desicji  can  be  stated  succinctly  in  one 
sentence;  an  intellii^ent,  choice  of  operating  band  is  sufficient.  The 
remainder  oi‘  this  section  jiresents  an  expanded  discussion  of  this  thoufdit. 

This  report  basically  constitutes  a  study  of  ttje  significance 
of  sipnal  spectral  cliaracteristics  upon  the  performance  of  a  conventional 
passive  detection  system.  Generally,  the  vicissitudes  of  the  sonar  environ¬ 
ment  are  such  that  a  detailed  knowledge  of  signal  characteristics  is  not 
possible.  Under  these  conditions,  an  Eckart  filter  implementation  is  not 
practically  feasible.  Tlie  results  in  this  report  suggest  two  possible  sub¬ 
optimum  systems:  the  first  is  a  conventional  pre-whitening  system  with 
a  bandwidtli  chosen  to  minimize  the  worst  case  loss  over  the  expected  raripv: 
of  signal  variations;  the  second  is  a  mismatched  (band-limited)  Eckart 
filter  with  a  transfer  function  again  chosen  to  minimize  worst  case  loss. 

Of  these  two,  the  first  is  probably  more  practical.  This  is  because  real¬ 
istic  estimates  of  signal  spectra  are  not  susceptible  to  a  simple  character¬ 
ization,  as  by  a  single  parameter  §.  Propagation  effects  can  markedly 
alter  the  behavior  of  it  would  be  difficult  to  make  a  plot  such 

as  Figure  H;  because  of  the  complicated  behavior  of  .Gi(f),  the  phrase 
"mismatched  Eckart  filter"  is  almost  meaningless.  It  is  likely  that  there 
would  exist  a  compromise  filter  that  provides  better  worst  case  performance 
than  the  pre-whitening  system,  Imt  t.he  likelihood  of  such  a  filter  being 
found  by  trial  and  error  methods  is  probably  small.-*  Further,  the  results 
in  this  report  imply  that  the  available  improvement  in  performance  is  rather 
small.  Consequently,  there  appears  to  be  considerable  justification  for 
being  content  wit.h  the  design  of  a  conventional  system  under  the  constraint 
that  the  shaping  filter  perform  only  the  function  of  pre-whitening;  the 

'<.;ee  Ap})endix  D(2)  . 
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Lhe  noise  and  band.limiiin(T;  the  input. 

Note  that  within  this  constraint  Wiere  are  still  considerable 
j’iains  to  be  achieved  by  makinr  intellip;ent  use  of  signal  spectral  infor- 
tration.  The  numerical  results  contained  in  this  report  (Figures  4  and  5) 
quite  possibly  are  not  applicable  to  the  realistic  sonar  situation,  but 
the  technique  of  analysis  is.  Thus,  one  could  apply  equations  1,  2,  and 
10  to  more  realistic  estimates  of  the  various  sonar  parameters  (T(f), 

Il(R,f),  N,(f),  etc.),  examining  the  influence  of  perturbations  in  these 
quantities  upon  system  perj'ormance  (the  outjiut  hNR)  with  different  operat¬ 
in'^  bands,  ouch  a  study  would  be  more  profitable  t.han  tfiat  carried  out 
here  since,  first,  the  spectra  would  be  more  representative  of  the  real 
world,  and  second ,  the  study  would  determine  both  ii  and  W,  i.e.,  ijoth  the 
location  and  widlli  of  the  operatiri!'-  hand,  f.ucli  a  comprehensive  study  may 
appear  formidable  and  len.  thy,  but  the  utilization  oi'  di'’-iLal  computer 
t.echniques  should  reduce  t,tie  effort  required  to  an  acccjit.ahle  level. 

Thus  we  return  once  more;  an  intelligent,  ciioice  of  operat.in!' 
band  is  sufficient.  ..tuc)i  a  conclusion  has  i]ndoub1,edly  been  bantered  about 
intuitively  for  some  l,imo,  so  that  some  would  question  tlie  value  of  the 
material  contained  in  this  report;  however,  the  quantitative  results  pre¬ 
sented  here  for  simple  spectra  should  provide  a  firmer  basis  for  such  an 
intuitive  judgement,  and  as  such,  provide  firmer  around  for  the  conventional 
system  designer  to  stand  on. 
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APPENDIX  A:  SNRo  EQUATION  FOR  SINGLE-SIDED  POWER  SPECTRA 
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The  derivation  of  equation  (9)  utilized  double-sided  power 
spectra;  however,  in  obtaining  numerical  results,  it  is  often  convenient 
to  work  with  single-sided  spectra.  The  necessary  modification  to  equa¬ 
tion  9  is  derived  below. 

The  double -sided  result  is  (Eqn.  9) 

Letting  a  super  A  denote  single-sided  spectra,  we  have  the  following 
conversions;  .  r 

3  sCflc/f 

“  o 

m  « ^  J 

I  o  /“o 

Thus 


since  N(f )  is  an  even  function.  Thus 

and  equation'  9  can  be  rewritten  for  single-sided  spectra  as 

=T  -  • 
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B«eau>e  of  Ite  siinpllclty,  e(tuation  (ll)  1b  particularly 
appealing.  It  clearly  reveals  the  small  signal  suppression  effect  — 
the  dependence  of  the  output  SNR  upon  the  square  of  the  Input  SNR  ••  • 
which  Is  characteristic  of  all  non-llnear  detectors  at  low  Input  SNR's 
(U,  pp  279“286;  lU,  pp  267  A  307).  It  also  Indicates  the  Importance  of 
the  "signal* processing"  parameters,  T  and  W.  The  general  result  In 
eqiuatlon  (9)  can  be  put  In  the  form  of  the  white  noise  result  In  equation 
(ll)  by  defining  an  appropriate  equivalent  noise  bandwidth,  viz. 

w  - 1  r 

Then  equation  9  can  be  written  as  ^e 

‘  f cS  P 


ritten  t 

Mr 


jf] 

irijcfjM 

j  =  7-Wt 


(61. fc) 

($/.  b) 


(BO 


where  S  and  N  denote  total  Input  signal  and  noise  powers  -  l.e.,  the 
result  of  Integrating  the  Input  power  spectra  S(f )  and  N(f). 

There  Is  precedence  for  Introducing  a  definition  such  as 
equation  Bl:  e.g.,  Blackman  and  Tukey  (16,  p  19  ff)  utilize  this  deflnl> 
tlon  of  equivalent  width  In  their  analysis  of  power  spectra  measurement. 
Although  W,  may  appear  to  be  on  obscure  method  of  defining  bandwidth, 
rewriting  and  B«nlpulatlng  equation  Bl  will  reveal  a  i^iysical  Interpre- 
tatlon  of  the  result.  Thus 
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Aside  from  the  factor  of  W,  is  the  ratio  of  the  area  under 
the  autocorrelation  (autoconvolution)  function  of  N(f)  to  the  value  of 
the  autocorrelation  function  evaluated  at  zero  shift.  Thus,  W,  can  be 
referred  to  as  the  equivalent  autocorrelation  width  of  N(f );  Bracewell 
(17;  PP  152-155)  presents  a  discussion  as  to  why  such  a  "width  mecsure" 
can  bo  advantageous  under  certain  conditions,  although  liis  measure  does 
not  Include  the  factor  of 

Tlie  concept  of  an  equivalent  autocorrelation  width  has  been 
freely  used  in  the  time  Jjm.-iin.  i.f-.  lor  F,  (x)  (9,  pp  2  <’.-24B) ;  there,  the 
width  is  referred  to  us  the  correlation  time,  t,  ^  for  the  process  x(t). 

The  actual  correlation  function,  R»(t),  is  replaced  by  a  rectangular 
function  of  width  t,  and  licight  K|((0).  It  can  readily  be  seen  that  an 
analogous  interpretation  is  valid  for  the  frequency  domain;  W«  is  one  of 
the  many  possible  —  though  perhaps  least  used  —  definitions  of  band- 
'-Idth. 

TJie  factor  of  \  my  be  botliersonie,  but  it  is  basically  Just  fall¬ 
out  from  the  raathemticc.  It  could,  uf  course,  be  arbitrarily  left  out 
of  the  definition  in  eqn.  (Hl),  and  a  factor  of  |  would  then  appear  in  eqn 
(R2).  It  is  more  convenient  and  intuitively  satisl'ying  to  leave  it  in 
the  definition  of  bandwldtli.  As  defined  in  (bl),  both  idetl  low  pass 
!ind  ideal  band  puss  (reuctungular  functions)  white  noise  of  bandwidth  W 
have  W,  =  W.  niraib'iTly,  a  triangular  spectrum  of  total  extent  W  has 
We  =  0.75  Otljer  examples  are  provided  in  both  Plackinan  and  Tukey 
(16)  and  Bracewell  (17). 
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A  convenient  way  to  determine  the  critical  value  of  §  is  to 
study  both  sides  of  the  extrema  equation,  equation  l8.  Define  the 
following  two  functions: 

-hi 

where  §  /  1,  p.  >  0.  Then  equation  l8  can  be  restated  as 


and  the  extrema  of  are  determined  by  the  intersections  of  f^(^i)  with 

It  is  readily  verified  by  substitution  that  U  =  0  is  always  a 
solution,  but  since  this  corresponds  to  zero  bandwidth,  it  is  not  the 
maximum  of  we  are  seektlng.  f^Ck*)  is  a  straight  line,  and 

if  5  <  it  has  a  negative  slope  with  an  intercept  of  1  at  u  "  0. 
Temporarily  restricting  ?  to  be  >  0,  g^{u)  is  a  "square  root"  type  of 
curve,  always  greater  than  one,  with  a  value  of  unity  at  ki  ■  0,  Consequent¬ 
ly,  f^(ki)  does  not  Intersect  g^(k‘)  at  any  value  of  k*  >  0.  Thus  infinite 
bandwidth  solutions  are  obtained  for  all  values  of  5  <  .i.  (The  negative 
§,  infinite  bandwidth  solutions  are  a  trivial  extension  of  the  above  argu¬ 
ments  . ) 

For  5  ■  ^,  f^(k*)  is  a  line  of  zero  slope,  at  a  height  of  1,  and 
g|(k*)  is  Vro.  Thus,  in  this  case,  there  is  no  intersection  (ignoring 
the  trivial  kk  ■  0  conincldence ),  and  the  infinite  bandwidth  solution  again 
results . 


With  5  >  f^(k*)  is  now  a  strai^t  line  with  positive  slope. 

With  ^  <  1,  g^(ki)  is  etlll  a  square  root  type  of  curve.  With  §  > 

f^(ki)  is  approaching  ®  as  k*— ♦  *  linearly  with  k^;  with  i  <  5  <  1,  g^(ki) 
is  approaching  ®  as  u— ^  ®  as  u^,  ^diich  is  less  than  linear  with  k*. 
Consequently,  f^(ki)  will  intercept  g^(ki)  before  k^  becomes  infinite. 

These  concepts  are  Indicated  graphically  in  the  accompanying 
figures  for  ?  •  0,  0.5,  0.7.  Trial  and  error  shows  that  the  solution 
for  5  ■  0.6  occurs  at  a  value  of  ki  approximately  equal  to  U5. 
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This  appendix  expands  and  clarifies  some  of  the  statements 
made  in  the  text.  The  statements  are  of  a  general  nature  and  rather 
lengthy,  and,  as  such  were  deemed  inappropriate  for  inclusion  in  the 
main  body.  The  numbers  of  the  headings  refer  to  the  footnote  references 
in  the  text. 

(l)  Tliat  Bayes  and  maximum  SNR  lead  to  the  same  optimum  con¬ 
figuration. 


This  statement  is  based  upon  a  number  of  approximations 
and  assumptions  which  may,  at  first  glanr-'',  appear  un¬ 
justifiable;  however,  experience  indicates  that  the 
statement,  though  perhaps  not  rigorous,  is  reasonably 
correct .  The  logic  leading  to  it  is  based  on  the 
following: 

a.  If  the  noise  into  the  square-law  device  is  bandlimited 
white  and  Gaussian,  then  the  equi-spaced  Nyquist  samples  of  a  segment 

T  seconds  long  are  independent  and  obey  the  chi-squared  distribution  at 
the  output  of  the  integrator  (summer  for  sampled  values).  (l,pp  77“78) 

b.  For  TW  >  50,  the  chi-square  distribution  can  be  satisfac¬ 
torily  approximated  as  Ouassian  (18).  This  is  simply  an  illustration 
of  the  Central  Limit  Theorem,  but  note  that  the  chi-square  distribution 
is  sufficiently  well  behaved  for  the  Gaussian  approximation  to  be  valid 
at  rather  small  values  of  TW. 


c.  If  the  output  of  the  processor  is  Gaussian,  the  statisti¬ 
cal  description  of  its  performance  is  readlLv  characterized  by  a  single 
parameter  d,  defined  as  d  • 

where  it  is  assumed  that  Oj  =  o  =  c.  (n)  =  Oz{3  +  II).  It  will  now  be  shown 
that  maximization  of  this  parameter  results  in  a  .Gayes  optimization  for 
Guassian  output  statistics. 

The  system  false  alarm  and  detection  probabilities  are 
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where  Zj  =  threshold  value  of  Z,  and  .  /a  ^  > 

f,  (2)  »  3“)  '  ey^  [.  13^^  j 

The  subscripts  0  and  1  refer  to  N  and  S  +  N,  respectively,  and  Zo  and  Zi 


are  the  output  means. 


Define  two  new  variables, 

^ 


where 


-  C^)"'  f  «' 

^•T 

'h  [**  e 

t  -  ?T.- 


-  c 


-  7  ^ 

-  X 


/  /  f 

“I 


'^IT 


If  the  PpA  is  specified  as  some  constant  value,  then  §ot  is  likewise  a 
constant,  say  Co,  so  that 


5ot  -- 


=  tC,  t  2, 


Pd  ^(v5?/f 
r 

'O’ 

Under  the  non-restrictive  assumption  Zj  >  Zo,  i.e.,  the  mean  with  S  +  N 
is  greater  than  the  mean  with  N,  the  probability  of  detection  is  maxi- 

ir.ized  by  maximizing  the  quantity 

d-  a— 

Maximizing  Pd  while  Pp*  is  lield  constant  is  Just  the  Neyman-Pearson 
criterion,  which  differs  from  the  Bayes  criterion  only  in  the  specifica¬ 
tion  of  the  threshold  of  tl.e  likelihood  ratio  (l,  pp  69-70).  Consequent¬ 
ly  maximizing  d,  or  equivalently,  maximizing  the  output  CNR,  satisfies 
tlic  Bayes  criterion  for  Ciaussion  output  statistics. 
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It  1b  important  to  realize  the  context  in  which  these  statements 
are  made.  In  order  for  the  maximization  of  d  to  have  meaning,  a  process¬ 
or  structure  must  be  specified;  one  cannot  derive  the  structure  from  d, 
but  can  optimally  choose  parameters  of  the  structure  by  maximizing  d. 

For  an  unstructured  problem  in  which  one  is  attempting  to  determine  the 
best  manner  to  process  raw  data,  resort  must  be  made  to  likelihood  ratio 
techniques  (the  Bayes  criterion). 

The  preceding  logic  is  then  applied  to  the  problem  at  hand  by 
stating  that  although  the  output  statistics  may  not  be  chi-square  (the 
noise  is  not  necessarily  white),  they  probably  possess  the  characteris¬ 
tic  of  the  chi-square  distribution  that  the  output  distribution  approaches 
Gaussian  for  reasonable  values  of  TW, ,  say  TWe  greater  than  100  or  so. 

Then,  by  sissiomlng  small  input  SNR's,  so  that  02(0  +  N)^a2(N),  one  can 
state  that  maximizing  d  (or  the  output  SNR)  satisfies  the  Bayes  criterion. 

(2)  That  the  best  one  can  do  where  highly  limited  a  priori  sig¬ 
nal  information  is  available  is  to  whiten  the  noise  and  bandllmit  the  input. 

The  work  contained  in  the  tody  of  this  report  can  be  considered  an 
an  attempt  to  Justify  this  statement.  The  discussion  in  Section  8  is 
directly  applicable;  as  stated  there,  it  is  likely  that  there  exists  a 
compromise  filter,  one  which  does  not  necessarily  whiten  the  noise  and 
would  have  a  spectral  shape  somehow  dependent  upon  the  expected  range  of 
signal  variations  --  analogous  to  the  mismatched  Echart  filter  for  sinple 
spectra  treated  in  Section  7  —  but  that  such  a  filter  would  be  difficult 
to  find  since  it  would  have  to  be  determined  by  trial  and  error  techniques. 
The  quantitative  results  in  this  report  indicate  that  the  gains  achievable 
from  such  a  filter  are  minimal,  and  that  a  bandllraited  pre-whitening  system 
should  prove  satisfactory. 

In  the  practiced  situation,  of  course,  one  is  never  going  to 
build  filters  that  actually  whiten  the  noise  or  sharply  bandlimlt  the  in¬ 
put.  The  spectral  fluctuations  of  the  input  noise  emd  considerations  of 
physical  realizability  theory  make  these  idealistic  concepts;  the  power 
of  the  dollar  over  system  realization  makes  their  approximations  doubtful. 
The  results  In  this  report  tend  to  substantiate  the  hope  that  these 
limitations  do  not  severely  compromise  performance. 
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